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Abstract 

Making use of an ansatz for the eigenfunction, we obtain exact closed form 
solutions to the Schrodinger equation with the inverse-power potential, V(r) = 
Ar~ 4 + Br~^ + Cr~ 2 + Dr^ 1 both in three dimensions and in two dimensions, 
where the parameters of the potential A, B, C, D satisfy some constraints. 
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1. Introduction 



The exact solutions to the fundamental dynamical equations play an important role 
in the different fields of physics. As far as the Schrodinger equation concerned, the exact 
solutions are possible only for the several potentials and some approximation meth- 
ods are frequently used to arrive at the solutions. The problem of the inverse-power 
potential, l/r n , has been widely carried out on the different fields of classic mechan- 
ics as well as on the quantum mechanics. For instance, the interatomic interaction 
potential in molecular physics [1-2], the inverse-power potentials V(r) = —Z 2 a/r A [3] 
(interaction between an ion and a neutral atom) and V(r) = —d^/r 3 [4] (interaction 
between a dipole di and another dipole d^) are often applied to explain the interaction 
between one matter and another one. The interaction in one-electron atoms, muonic 
and hadronic and Rydberg atoms also requires considering the inverse-power potentials 
[5]. Indeed, the interaction potentials mentioned above are only special cases of the 
inverse-power potential when some parameters of the potential vanishes. 

The reason we write this paper is as follows. On the one hand, Ozcelik and Simsek 
discussed this potential in the three-dimensional spaces [6]. They obtained the eigen- 
values and eigenfunctions for the arbitrary node. Simultaneously, the corresponding 
constraints on the parameters of the potential were obtained. Unfortunately, they did 
not find that it is impossible to discuss the higher order excited state except for the 
ground state. In the later discussion, we will draw this conclusion and find some essen- 
tial mistakes occurred in their calculations even for the ground state. We recalculate 
the solutions to the Schrodinger equation with this potential in three dimensions fol- 
lowing their idea and correct their mistakes. On the other hand, with the advent of 
growth technique for the realization of the semiconductor quantum wells, the quantum 
mechanics of low-dimensional systems has become a major research field. Almost all 
of the computational techniques developed for the three-dimensional problems have al- 
ready been extended to lower dimensions. Therefore, we generalize this method to the 
two-dimensional Schrodinger equation because of the wide interest in lower-dimensional 
fields theory. Besides, we has succeeded in dealing with the Schodinger equation with 
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the anharmonic potentials, such as singular potential both in two dimensions and in 
three dimensions [9, 10], the sextic potential [11], the octic potential [12] and the Mie- 
type potential [13] by this method. We now attempt to study the Schodinger equation 
with the inverse-power potential by the same way both in three dimensions and in two 
dimensions. 

This paper is organized as follows. In section 2, we study the three-dimensional 
Schrodinger equation with this potential using an ansatz for the eigenfunctions. The 
study of the two-dimensional Schrodinger equation with this potential will be discussed 
in section 3. The figures for the unnormalized radial functions are plotted in the last 
section. 

2. Solutions in three dimensions 

Throughout this paper the natural unit H — 1 and /i — 1/2 are employed. Consider 
the Schrodinger equation 

-V 2 ^ + V(r)ip = Elf), (1) 
where here and hereafter the potential 

V{r) = Ar- 4 + Br~ 3 + Cr~ 2 + Dr' 1 , A > 0, D < 0. (2) 

Let 

^(r,e, V )=r- 1 R e (r)Y em (0,v), (3) 

where i and E denote the angular momentum and the energy, respectively, and the 
radial wave function Ri{r) satisfies 



d 2 Ri(r) + 



Mr) = 0. (4) 



dr 2 

Ozcelik and Simsek [6] make an ansatz for the ground state 

BP l (r)=exp\g(r)], (5) 

where 

g (r) = - + br + clnr, a < 0, b < 0. (6) 

r 
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After calculating, one can obtain the following equation 



d 2 R°e(r) 
dr 2 



d 2 g(r) I dg{r 
dr 2 



+ 



dr 



R°Ar) 



(7) 



Compare Eq. (7) with Eq. (4) and obtain the following sets of equations 



a 2 = A, b 2 = -E, 



2bc = D, 2a(l - c) = B, 



C + 



+ 1) -4 



26a — c. 



(8a) 

(86) 
(8c) 



It is not difficult to obtain the value of the parameter a from Eq. (8a) written as 
a = ±y/A. In order to retain the well-behaved solution at r — > and at r — > oo, they 
choose negative sign in a, i. e. a = —\f~A. According to this choice, they arrive at a 
constraint on the parameters of the potential from Eq. (8c) written as 



B 



+ 



2AD 



£(£+1). 



(9) 



4A 2y/A B + 
Then the energy is read as 

E% = [C + £(£ + 1) ± tJ[C + £(£+1)] 2 -2Bd} 2 . (10) 

It is readily to find that Eq. (10) is a wrong result. From Eqs. (6) and (8b), as we 
know, since the parameter b is negative, when we calculate the energy E from Eq. (8a), 
we only take the b as a negative value, so that Eq. (10) only takes the negative sign. 
Actually, it is not difficult to obtain the corresponding values of the parameters for the 
g(r) from Eq. (8), i. e. 

B + 2^A , D^fA 



2^/A ' " B + 2^A' 
The eigenvalue E, however, will be simply expressed as from Eq. (8a) 

E _ AD 2 

B 2 + AA + AB\/~A' 

The corresponding eigenfunction Eq. (5) can now be read as 



11) 



(12) 



R° e = N r c exp 



-a + br 

r 



(13) 



where N is the normalized constant and here and hereafter the parameters a, 6 and c 
are given above. 

After their discussing the ground state, Ozcelik and Simsek continue to study the 
first excited state. They make the ansatz for the first excited state, 

Rl(r) = f(r)eM9(r)], (14) 

where g(r) is the same as Eq. (6) and f(r) = r — a±, where a\ is a constant. For 
short, it is readily to find from Eq. (14) that the radial wave function R\(r) satisfies 
the following equation 



g(r)" + (g(r)'r + 



,. 2 , , f(r)" + 2g(ryf(ry 



R]{r) = 0, (15) 



f(r) 

where the prime denotes the derivative of the radial wave function with respect to the 
variable r. 

Compare Eq. (15) with Eq. (4) and obtain the following sets of equations 

—26 — 26c + D + b 2 ai + ea x = 0, (16a) 

-b 2 -E = 0, a 2 ai - Aa x = 0, (166) 

—a 2 + A + 2act\ — Bol\ — 2acct\ = 0, (16c) 

2a6 - c - c 2 + C + i{i + 1) + 26c«i - Da x = 0, (16d) 

B + 2ac — 2aboi\ — ca.\ + c 2 o;i — Ca\ — ta.\ — £ 2 ai = 0, (16e) 
it is not hard to obtain the following sets of equations from Eqs. (16a-16c) 

E = -b 2 , a 2 = A, c= B + 2 / ^ J , (17a) 

2VA V ; 

t= ^ (176) 

where the constant a\ ^ and it is determined by Eqs. (16d) and (16e). Furthermore, 
it is evident to find that Eq. (17b) does not coincide with Eq. (11) with respect to the 
same parameter 6, which will lead to the their wrong calculation for the first excited 
state. In fact, they obtained two different relations during their calculation through 



the compared equation, i. e. D — 2bc (see Eq. (9) in [6]) and D = 2b(c + 1) (see Eq. 
(16) in [6]). The parameter D does not exist if the parameter b is not equal to zero. It 
is another main mistaken that arises their wrong result, that's to say, it is impossible 
to discuss the first excited state for the Schrodinger equation by this method. We only 
discuss the ground state by this simpler ansatz method as mentioned above. 

As a matter of fact, the normalized constants No can be calculated in principle from 
the normalized relation 

/ \R° e \ 2 dr = 1. (18) 
Jo 

In the course of calculation, making use of the standard integral [14](ReAi > 0, ReA 2 > 
and Rez/ > 0) 

K v (2y/X^\ 2 ), (19) 

which implies 

N ° = _2(f) 2 ^ i K 2c+1 (4 v / ^) 
where the values of the parameters b, c and a are given by Eq. (11) and — v^4, respec- 
tively. The figure 1 for the unnormalized radial eigenfunction in three dimensions is 
plotted in the last section. 

3. Solutions in tow dimensions 

We now generalize this method to the two-dimensional Schrodinger equation. Con- 
sider Schrodinger equation with a potential V(r) that depends only on the distance r 
from the origin 

^ = {ll r lr + * + V{r)i) = E *' (21) 

Let 

$(r, ip) = r-^XWe^, m = 0, 1, 2, . . . , (22) 
where the radial wave function R m (r) satisfies the following radial equation 

R m (r) = 0, (23) 

6 



(20) 



d 2 R m (r) 
dr 2 



+ 



E - V(r) 



m 



1/4 



where m and E denote the angular momentum and energy, respectively. For the 
solution of Eq. (23), we make an ansatz [6-13] for the ground state 



where 



R° m (r) = exp[g m {r)], 



9m{r) = — + b l r + c l lnr. 

r 



(24) 



(25) 



After calculating, we arrive at the following equation 



dr 2 



d 2 g m {r) , / dg m (rY 



dr 2 



+ 



dr 



R° m (r) = 0. 



(26) 



Compare Eq. (26) with Eq. (23) and obtain the following sets of equations 

a\ = A, b\ = —E, 

2b lCl = D, 2ai(l - ci) = B, 

C + rr? — - = c? — 2b\a\ — c\. 
4 



(27a) 

(276) 
(27c) 



It is not difficult to obtain the values of the parameters a\ from Eq. (27a) written 
as di = ±v^4. Likely, in order to retain the well-behaved solution at r — > and at 
r — > oo, we choose negative sign in ai, i. e. ai = —y/A. According to this choice, Eq. 
(27b) will give the other parameter values as 

B + 2\/A , DVA 



Cl = 



2VA B + 2VA 

Besides, it is readily to obtain from Eq. (27c) that 



(28) 



„ B 2 B 2AD 
C = — + — = + 



(m 2 -l/4), 



(29) 



4A 2y/A B + 2y/A 

which is the constraint on the parameters for the two-dimensional Schrodinger equation 
with the inverse-power potential. 

The eigenvalue E, however, will be given by Eq. (27a) as 

AD 2 



E 



B 2 + AA + AB\/A' 



(30) 
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The corresponding eigenfunction Eq. (24) can now be read as 




Nr Cl exp -a\ + b±r , 



r 



(31) 



Similarly, the normalized constants iV can be calculated in principle from the nor- 
malized relation 



where the values of the parameters a 1; b\ and C\ are given above. The figure 2 for the 
unnormalized radial eigenfunction in two dimensions is plotted in the last section. 

Considering the values of the parameters of the potential, we fix them as follows. 
The values of parameters A, C, D are first fixed, for example A = 4.0, C = 2.0 and 
D = —2.0, the value of the parameter B is given by Eq. (10) and Eq. (29) for the 
cases both in three dimensions and two dimensions for £ = and m — 0, respectively. 
By this way, the parameter B turns out to B = 5.87 in three dimensions and B = 5.65 
in two dimensions, respectively. The ground state energy corresponding to these values 
are obtained as E = —0.164 for the case in three dimensions and E = —0.172 for the 
case in two dimensions. Actually, when we study the properties of the ground state, 
as we know, the unnormalized radial wave functions do not affect the main features of 
the wave functions. We have plotted the unnormalized radial wave functions in figures 
1 and 2 for the cases both in three dimensions and in two dimensions, respectively. 
With respect to figures 1 and 2, it is easy to find that they are similar to each other, 
which stems from the same values of the angular momentum I = and m — 0. They 
will be different if we take the different values of the angular momentum in the course 
of calculations. 

In conclusion, we obtain the exact analytic solutions to the Schrodinger equation 
with the inverse-power potential V(r) = Ar~ 4 + Br~ 3 + Cr~ 2 + Dr^ 1 using a simpler 
ansatz for the eigenfunction both in three dimensions and in two dimensions, and 




(32) 



According to Eq. (19), we can obtain 



N = 



1 



(33) 



2(fJ-)^ K 2ci+l {A^h) 
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simultaneously the constrains on the parameters of the potential are arrived at from 
the compared equations. Finally, we remark that this simple and intuitive method 
can be generalized to other potential. The study of the Schrodinger equation with the 
asymmetric potential is in progress. 
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5 10 15 20 25 

Figure 1: The ground state wave functions in three dimensions as a function of r for 
the potential (2) with the values A = 4.0, B = 5.87, C = 2.0 and D = -2.0. The y-axis 
denotes the values of wave functions and the rr-axis denotes the variable r. 
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Figure 2: The ground state wave functions in two dimensions as a function of r for the 
potential (2) with the values A = 4.0, B = 5.65, C = 2.0 and D = -2.0. The y-axis 
denotes the values of wave functions and the x-axis denotes the variable r. 
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